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METHOD FOR MEASURING PHYSICAL PARAMETERS OF AT LEAST 
ONE MICROMETRIC OR NANOMETRIC DIMENSIONAL PHASE IN A 

COMPOSITE SYSTEM 

The subject of the present invention is a method of 
5 determining at least one mechanical parameter of at 
least one material in a composition system comprising 
at least two distinct phases. The method applies more 
particularly to the case in which the phase to be 
characterized is of a microscale or nanoscale 
10 structure, and at least one other material constitutes 
a substrate or a matrix. For example, the method 
applies to the characterization of a thin film 
deposited on a substrate or of inclusions, filaments or 
fibers in a matrix. 

15 

The expression "mechanical parameters" is used here to 
indicate both mechanical and thermomechanical 
properties, such as the Young's modulus, the Poisson's 
ratio, the thermal expansion coefficient, etc., and 
20 strain and/or stress states. 

Knowing the mechanical properties of a system or 
physical device allows its operation to be optimized. 
Any material, whatever it is, is subjected to external 

25 stresses and is necessary to be able to determine its 
resistance to such stresses. It is therefore important 
to completely determine the elastic properties of a 
system and in particular those of one or more layers 
placed on the surface of a substrate. This is because 

30 the properties of this layer (or these layers 
considered as a subsystem) differ appreciably from the 
properties of the films taken individually. 

The elastic properties of materials are involved in 
35 many fields of application: coating of mechanical 
parts, structural deformation, etc. They play an 
increasingly important role in the fabrication of 
electronic circuits: since the size of devices is 



decreasing, the strains generated at the interface of 
the various parts of the device are correspondingly 
higher . 

The methods already developed for measuring the elastic 
constants and expansion coefficients are generally 
carried out on a macroscale. 

The conventional methods are: 

A) deflection or deformation of macroscopic or 
smaller specimens, see for example: " Measurement of 
elastic modulus,, Poisson ratio and coefficient of 
thermal expansion of on-wafer submicron films" by Jie- 
Hua Zhao, Todd Ryan, Paul S. Ho, Andrew J. McKerrow and 
Wei-Yan Shih, PRC/MER2 . 206, Journal of Applied Physics 

(1999) , 85/(9) , 6421-6424; or "Elastic modulus 
measurements by three methods on a 2-D laminated 
continuous fiber-reinforced ceramic matrix composite" 
by S.T. Gonczy and M.G. Jenkins, Ceramic Transactions 

(2000) , 103 (Advances in Ceramic-Matrix Composites V) , 
541-547; or else "Analysis of residual stress in cubic 
boron nitride thin films using micromachined cantilever 
beams" by G.F. Cardinale, D. G . Howitt, K.F. McCarty, 
D.L. Medlin, P.B. Mirkarimi and N.R. Moody, Diamond and 
Related Materials (1996), 5(11) , 1295-1302; 

B) electroacoustic waves: for example, "Hardness 
and Young's modulus of high-quality cubic boron nitride 
films grown by chemical vapor deposition" by X. Jiang, 
J. Philip, W.J. Zhang, P. Hess and S. Matsumoto, 
Journal of Applied Physics (2003), 93 (3) , 1515-1519 or 
" Thin- film elastic-property measurements with laser- 
ultrasonic SAW spectrometry" by D.C. Hurley, 
V.K. Tewary and A.J. Richards, Thin Solid Films (2001), 
398-399) , 326-330; and 

C) nanoindentation : for example, "Determination 
of mechanical film properties of a bilayer system due 
to elastic indentation measurements with a spherical 
indenter" by T. Chudoba, N. Schwarzer, F. Richter and 
U. Beck in Thin Solid Films (2000), 377-378 , 366-372 or 



"Hardness and elastic modulus measurements of A1N and 
TiN sub-micron thin films using the continuous 
stiffness measurement technique with FEM analysis" by 
T.A. Rawdanowicz, J. Sankar, J. Narayan and V. Godbole, 
Materials Research Society Symposium Proceedings 
(2000), 594 (Thin Films-Stresses and Mechanical 
Properties VIII), 507-512. 



There are also more recent methods for measuring the 
mechanical properties of small objects, for example 
atomic force microscopy (AFM) described for example in 
"Measurements of elastic properties of ultra-thin 
diamond-like carbon coatings using atomic force 
acoustic microscopy" by S. Amelio, A.V. Goldade, 
U. Rabe, V. Scherer, B. Bhushan and W. Arnold, Thin 
Solid Films (2001), 392 (1) , 75-84. 

The object of the invention is to measure one or more 
mechanical parameters, especially the elastic constants 
and the thermal expansion coefficients of a material 
and the strains that this material generates when it is 
combined with other materials. Its specificity is to be 
able to measure such parameters under conditions close 
to their application conditions: thin or thick films, 
inhomogeneous layers, discontinuous "layers" or 
precipitates/inclusions, fibers or filaments, and boxes 
(the term more specific to semiconductor materials in 
which charges are localized in these boxes formed from 
a second material different from the substrate) . 
Current knowledge prevents one from asserting that 
microscopically measured parameters can be used in 
structures of small dimensions, that is to say on a 
microscale or nanoscale. 



The basic idea of the invention consists in using a 
phenomenon that is well known but is generally 
considered as a drawback, which it is sought to avoid 
or neglect: in a specimen having at least one 
sufficiently small dimension, for example on the micron 
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or submicron scale, owing to the proximity of the free 
surface, the strains built up in the initial device are 
relaxed (see the aforementioned documents by 
D. Perovic, G. Weatherly and D. Houghton, Phil. Mag. 
5 64, (1991)). According to the invention, it is 
unnecessary to make direct measurements on the material 
to be characterized (something that is not always 
possible or easy) : by measuring the deformation 
parameters of the substrate or matrix, the origin of 

10 which is the relaxation of the strains in the specimen, 
it is possible to obtain the parameters of the phase 
with the microscale or nanoscale structure. Any 
measurements made on this phase, when they are 
technically possible, provide additional information 

15 but they are not in general essential. 

The expression "deformation parameter" is understood to 
mean both pure deformations, usually indicated by the 
symbol e and which, from the microscopic standpoint, 
20 correspond to changes in the crystal parameters, and 
local rotations of crystal axes, indicated by the 
symbol p. 

The invention therefore relates to a method for 
25 determining at least one mechanical parameter of at 
least one material in a composite system comprising at 
least two distinct phases, characterized in that it 
comprises : 

a) the production of at least one specimen 
30 comprising a first part of a first phase and a second 
part of a second phase, the second part consisting of 
the material to be characterized, the specimen having 
at least one dimension small enough to allow the 
strains in said specimen to be relaxed; 
35 b) the measurement, on said specimen, of at 

least one deformation parameter of at least said first 
phase, in correspondence with a plurality of points 
lying at different distances from an interface between 
said first and second phases; and 



c) the 
deformation 
parameter of 



determination, 
parameter, of at 
said second phase. 



from at least said 
least one mechanical 



According to particular embodiments: 
said method comprises: 

i) the production of a plurality of specimens 
that differ from one another in respect of at least one 
geometrical property; 

ii) the implementation of step b) on each of said 
specimens; and 

iii) the use in step c) of the measurements made 
on said plurality of specimens; 

step b is repeated at at least two different 
temperatures for at least one of said specimens; 
step c) comprises: 

i) the modeling of the strain relaxation in said 
specimen (L) using a first estimate of at least one 
mechanical property of the material of said second 
phase (B) ; 

ii) the comparison of the measurement results of 
step b) with those of said modeling; and 

iii) the modification of said estimate of at least 
one mechanical property of the material of said second 
phase and the reiteration of substeps i) to iii) until 
the difference between said measurement results and the 
modeling results is minimized. 

Within this context, the term "modeling" covers both 
numerical simulations and approximate analytical 
models. Advantageously, a finite-element numerical 
simulation may be used; 

said composite system is chosen from among: a 
substrate having a continuous layer on its surface; a 
substrate having metallization bands or islands on its 
surface; a layer with a zone included in the substrate; 
a transistor; a layer on the inside of a substrate; a 
matrix containing inclusions; fibers or filaments; 

- said specimen has at least one microscale or 
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nano scale dimension; 

said specimen is a lamella having two 
approximately parallel faces lying approximately 
perpendicular to the interface between said first and 
5 second phases, in which case step b) is advantageously 
repeated for a plurality of lamella (e) of different 
thicknesses; 

according to one alternative embodiment, said 
specimen is a lamella placed at an angle to the 
10 interface between said first and second phases, in 
which case step b) is advantageously repeated for a 
plurality of lamella (e) lying at different angles to 
the interface between said first and second phases; 

according to another alternative embodiment, 
15 said specimen is a wedge-shaped lamella having two 
faces making an angle between them, in which case step 
b) is advantageously repeated for a plurality of 
lamella (e) having two faces making different angles 
between them; 

20 - the measurements provided in step b) are 

carried out by diffraction of a convergent electron 
beam; 

step b) includes the observation of Holz 
lines for at least one crystallographic plane of said 
25 first phase and the determination of at least one 
parameter from among: the width of said Holz lines, 
their position and their internal structure; 

step b) comprises the determination of at 
least the width of at least some of said Holz lines and 
30 the calculation, for each of them, of a maximum 
rotation p max along the axis of the electron beam; and 

step c) involves the plotting of at least one 
curve representing a said maximum rotation as a 
function of the distance relative to the interface 
35 between said first and second phases. 

According to a preferred embodiment, step c) also 
involves, by simulation, the plotting of curves 
representing the maximum rotation p max as a function of 



the distance relative to the interface between said 
first and second phases for possible values of Young's 
modulus and/or Poisson's ratio of the material of said 
second phase, and also the minimization of the 
difference between the simulated curves and the 
experimental curves in order to determine the Young 1 s 
modulus and/or the Poisson's ratio of the material of 
said second phase. 

The invention will be more clearly understood on 
reading the following description in conjunction with 
the drawings in which: 

figures la to If illustrate various 
configurations of the device to be studied: figure la, 
a single film; figure lb, metallization bands or 
islands; figure lc, a film with a zone included in the 
substrate; figure Id, a transistor; figure le, a thin 
film inside a substrate; and figure If, fibers in a 
matrix. In general, A denotes the substrate or matrix 
and B denotes the microscale or nanoscale phase to be 
characterized; 

figures 2a and 2b show a sectional view of a 
transistor having a substrate made of silicon Si, which 
there is an NiSi film with a thickness e = 20 nm and an 
Si 3 N 4 film. A tungsten contact constitutes the drain 
contact D. The zone analyzed (figure 2b) is indicated 
by the arrow F. The thin lamella shown in figure 2b has 
a thickness t (which can be varied) . The mean incident 
beam of electrons for the CBED patterns (direction z 0 ) 
is taken along the y 2 axis that makes an angle y with 
the yi axis normal to the lamella. The directions and 
the corresponding crystallographic axes are illustrated 
in figure 2b. Note, CBED stands for Convergent Beam 
Electron Diffraction; 

figure 3a is a montage of photographs 
illustrating five CBED diffraction patterns chosen from 
around fifty images actually claimed, along a straight 
line perpendicular to the surface (the -z 2 direction) 
at a distance of 155 nm from the drain contact D. For 
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each pattern, the distance -z 2 (in nanometers) from the 
surface of the substrate (the interface with the thin 
film) and the measured angle A0 are indicated, the 
curve in figure 3b showing the variation in the 
5 calculated angle A9 as a function of -z 2 ; 

figure 4a shows a CBED diffraction pattern 
taken along a [230] zone axis in unstrained silicon 
with a large depth z 2 = 300 nm (x 2 // [3, -2, 0] , 
y 2 //[2,3,0], z 2 //[001]; 

10 - figure 4b is the pattern 4a in which Holz lines 

positioned by the "JEMS" software have been adjusted. 
The position of these lines makes it possible to 
determine, accurately, the semiconvergence angle of the 
incident beam and the accelerating voltage (the method 

15 . here is similar to that used in the "STREAM" project) ; 

figure 5a is a CBED pattern along the [230] 
zone axis in deformed silicon (CBED pattern of the C 
type taken at the depth z 2 = 140 nm) ; 

in figure 5b, lines obtained from three 

20 simulations with the "JEMS" software have been 
superimposed on this pattern: 

• the first system of lines represents the 
simulation of the diffraction pattern of a perfect 
silicon crystal tilted through by +Q ma x relative to the 

25 x 2 //[3-20] axis, 

• the second system of lines represents the 
simulation of the diffraction pattern of a perfect 
silicon crystal tilted through by -9 max relative to the 
x 2 //[3-20] axis, and 

30 • the third system of lines (dashed) represents 

the simulation of the diffraction pattern of the 
perfect silicon crystal but not tilted. The angle 20 max 
is denoted by AG; 

figures 6a to 6d illustrate the feasibility of 

35 minimizing the elastic coefficients of the material. 
Three parameters have been minimized manually, by trial 
and error, namely the coherence temperature T 0 (an 
imaginary temperature at which the two materials would 
be coherent), the Young's modulus E and the Poisson's 
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ratios v. 

The squares represent the experimental measurements of 
the angle AO. The curves with many points are obtained 
5 from finite-element simulations and represent the angle 
0.98AP, where AG = 2p max . The minimization criterion % 
determines the difference between the measured angles 
AG and the calculated angles 0.98AP, for various depths 
z 2 below this silicide and for different lamella 
10 thicknesses t; 

figure 6a: lamella thickness t = 300 nm, 
initial values of the parameters taken from the 
literature: T = 410°C; E = 150 GPa; v = 0.1; 

figure 6b: lamella thickness t = 300 nm, 
15 initial values of the parameters taken from the 

literature: T 0 = 410°C; E = 150 GPA; v = 0.1; 

figure 6c: lamella thickness t = 320 nm, 
final values after partial manual minimization: 
T 0 = 430°C; E - 115 GPa; v = 0.288; 
20 - figure 6d: lamella thickness t = 320 nm, 

final values after partial manual minimization: 
T 0 = 430°C; E = 115 GPa; v = 0.288; 

figure 7a is a CBED diffraction pattern 
obtained for a lamella thickness t = 320 nm at the 
25 distance -z 2 = 84 nm from the surface; 

figure 7b is a profile obtained over the (537) 
Holz band of figure 7a. The angular difference AG g is 
marked by two vertical lines and represents the width 
of the Holz band; 
30 - figure 7c is a simulation of the profile of 

figure 7b, which shows that it is possible to reproduce 
the broadening of the Holz lines and the variations in 
intensity 1 - A* (G) in the Holz bands. This simulation 

used the results of the finite-element calculation, 
35 which is an illustration of point iv (optimized elastic 
constants and displacement R(yiZi)). AG g is equal to 
0.11°, which corresponds to rotation of the x 2 axis 
through an angle AG = 0.14°. In the preliminary 
calculations, the widths AG g of the x 2 //[320] axis were 



ignored. A more refined minimization algorithm would 
seek to reproduce the curve similar to 7b in their 
entirety, with their oscillations and not to reproduce 
merely their width A0 g over the angle AG; 

figures 8a to 8d illustrate a second 
experimental example, in which the silicon substrate A 
is surmounted by an Si (1 _ x) Ge x film then by an Si film on 
the surface. Figure 8a shows the thin lamella of 
thickness t. Figures 8b to 8d are CBED diffraction 
patterns in the [230] direction within the Si(i- x) Ge x 
film, in the Si substrate far from the deformed zone, 
and in the Si substrate close to the deformed zone, 
respectively; and 

figure 9 shows a flow chart for one way of 
implementing the method of the invention. 

The detailed description that follows relates to the 
particular case in which the phase to be characterized 
is a thin film, the specimens used are in the form of a 
thin lamella of different thicknesses, the deformation 
parameter in question is the rotation of the crystal 
axes in a plane perpendicular to the film, and the 
measurement technique adopted is CBED (convergent beam 
electron diffraction) . 

In this case, the method of the invention comprises: 

a) the production of a lamella L with a 
sufficiently small thickness t and having two 
approximately parallel faces being approximately 
perpendicular to said substrate surface; 

b) the measurement, on said lamella, of at least 
one deformation parameter of the substrate at the 
various depths from the surface; deformation parameter 
is also understood to mean rotation parameter of the 
substrate; and 

c) the determination of at least one mechanical 
parameter of said film from at least said 
deformation/rotation parameter . 



The method may include the production of several 
lamella (e) of different thicknesses and the 
implementation of step b) on each of said lamella (e). 
For at least one said lamella, step b) may be repeated 
at at least two different temperatures. 

Advantageously, said measurement is carried out by 
generating, for points on the substrate located at 
different depths, CBED diffraction patterns for a 
convergent electron beam of axis Z 0 disoriented 
relative to the normal to said lamella, said patterns 
having Holz lines or bands. The determination at c) may 
then include the measurement of the width of the Holz 
lines of at least some of said patterns, for at least 
one crystallographic plane of the substrate. It is 
possible to calculate, from the width of these Holz 
lines, a maximum rotation p max along the electron beam 
axis for each pattern. This rotation is induced by the 
layer (or layers) placed on the substrate and said 
angle characterizes the properties thereof. 

Next, a curve can be plotted that represents said 
rotation p max as a function of the depth at which said 
patterns were obtained. It is then possible to plot, by 
simulation, at least one curve representing this 
rotation as a function of the depth for possible values 
of the parameters to be extracted, for example Young 1 s 
modulus and/or Poisson's ratio of the film to be 
analyzed in the case of an isotropic approximation and 
to minimize the difference between at least one 
simulated curve and a corresponding experimental curve, 
in order to determine these parameters. A similar 
technique may be used in the case of an isotropic 
modeling with parameters known to those skilled in the 
art. Values of the parameters that resulted in the 
simulated curve closest to the experimental curve are 
taken . 

The aforementioned four points of the invention will be 
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presented below. 

point ( i ) : the thin lamella of controlled geometry 
is extracted from or thinned in the device. A parallel- 
5 faced lamella is preferable but not absolutely 
essential. A slight angle may be present. A focused ion 
beam ( FIB) was used, but alternative methods, that are 
conventional in specimen preparation for electron 
microscopy, may be used (mechanical thinning, cleavage, 
10 etc.), but the FIB technique has the advantages of 
being rapid, of not mechanically disturbing the system 
or device, and of fully controlling the operations. 

point ( ii ) : to measure the deformations, a 

15 convergent electron beam (CBED) is used for example in 
scanning mode. A nanoscale electron beam is focused 
onto various points of the specimen so as to completely 
map the deformations of the thin lamella (figure 3) . A 
CBED diffraction pattern is therefore obtained at each 

20 point. This pattern comprises many Holz lines (more 
than 10 lines) , each corresponding to a 
crystallographic plane indexed by a vector g of the 
reciprocal lattice. It is preferable to acquire the 
CBED patterns in the substrate (part A) , since part B 

25 is generally too thin or too defective for being able 
to acquire good diffraction patterns in this zone b. 
However, for certain well-crystallized systems, 
patterns in zones a and b may be obtained. This would 
for example be the case if the layer a is a silicon 

30 layer and film b is composed of two thin films, namely 
an Sii_ x Ge x film (with a low Ge composition x, for 
example 10%) and an Si film. 

For a given lamella thickness and a given temperature 
35 of observation, an entire series of CBED patterns is 
produced so as to determine the crystal parameters and 
rotations of the substrate. 



point ( iii ) : by varying the lamella thickness (the 
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FIB technique allows this to be done easily) and the 
temperature of observation in the microscope, a whole 
series of experimental data points is obtained which 
make it possible (point iv) to calculate the elastic 
5 constants, the expansion coefficient of part B and the 
strains in the unthinned initial device. If part B is 
composed of a homogeneous film of a given material, it 
is unnecessary for the measurements to be carried out 
on the same lamella thinned in succession to different 

10 thicknesses. However, working on a single lamella the 
precision is increased, and this is essential if the 
system consists of a single nanosystem (for example a 
transistor) . A single lamella thickness does not allow 
all the constants of the material to be correctly 

15 calculated. A single lamella thickness gives 
information only about the strains in the material 
(this is the case of the study of the curvature of 
semiconductor wafers via the Stoney formula) . The 
strains are partly relaxed by a curvature of the 

20 substrate (see for example: "Measurement of elastic 
modulus, Poisson ratio, and coefficient of thermal 
expansion of on-wafer submicron films" by Jie-Hua Zhao, 
Todd Ryan, Paul S. Ho., Andrew J. McKerrow and Wei- 
Yan Shih, Journal of Applied Physics (1999), 85(9) , 

25 6421-6424) . 

point (iv) : simulations using the theory of 
elasticity are then performed in order to reproduce the 
experimental results. In complex systems, only finite- 
30 element calculations may be carried out. In simpler 
systems (part B = thin film) , an analytical formula may 
be implemented. The simulations reproduce the following 
phenomena : 

in the initial device, the various parts are 
35 strained or partially strained - the fact of extracting 
a thin lamella from the initial device relaxes the 
strains in the form of a rotation/change of crystal 
parameters and it is the observation and simulation of 
this strain relaxation that allows us to determine the 
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parameters of part B. 

The invention will be widely applicable in the surface 
treatment of mechanical parts, optimization of 
5 electronic circuits (metal contact, oxide film, etc.) 
or devices where the presence of two different 
materials necessarily creates mechanical strains. 

The method according to the invention is novel although 

10 it does make use of techniques or physical effects that 
are well known: 

the FIB (Focused Ion Beam) technique for 
specimen preparation; 

the convergent beam technique (CBED: Conversion 

15 Beam Electron Diffraction) , a particular technique used 
in electron microscopy. The STREAM project (see for 
example the publication "Software for Automation of 
TEM/CBED Methodology for strain determination" , IST- 
1999-10341 STREAM Consortium - Deliverable D23) uses a 

20 similar technique (FIB, convergent beam, simulation), 
but the technique essentially measures variations in 
the crystal parameters, whereas the method according to 
the invention is interested mainly in the local 
rotations of the crystal lattice. In addition, the aim 

25 of the STREAM project is not to measure elastic 
constants but to measure containted deformations 
strains in integrated circuits. 

The STREAM project measured only changes in crystal 
30 parameters far from the two parts A and B of the device 
and neglected strain relaxation in the thin lamella. 
According to the invention, detecting the rotation of 
the crystal lattice makes the method more rapid and 
more precise, and allows the interface between parts A 
35 and B to be approached. 

The power of the method according to the invention has 
been demonstrated by analyzing the strains produced by 
an NiSi film in an integrated circuit (specimen 1 in 
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figure 2a) . This system is relatively complex, as there 

are several electrical contacts and several materials 

present. For a rigorous treatment (which is quite 

possible), all the components of the system would have 
5 to be taken into account. 

In a first analysis, it has been assumed that part B 
can be considered as a thin film of infinite lateral 
extent surmounted by the atmosphere (the layers 
10 situated above the NiSi film were experimentally 
removed, and the thickness of the thin film was 
measured to be equal to t = 20 nm (see figure 2b) . 

The precise way in which the four points of part 3 are 
15 carried out will be described below. 

To do this, it is necessary to define various geometric 
coordinate systems . 

20 R 0 = (x 0 ,yo,zo) denotes the geometric coordinate system 
tied to the microscope. The z 0 axis is parallel to the 
optical axis of the microscope, defined as being the 
mean direction along which the electrons propagate 
before the specimen. The images or the diffraction 

25 patterns are recorded in the (x 0 ,yo) plane. 

Various coordinate systems R c = (x c ,y c ,z c ), 

Ri = (xi,yi,zi) and R 2 = (x 2 ,y2/Z 2 ), tied to the crystal 
structure of the silicon substrate (part A) , are also 
30 defined, where: 

Xc=nO0j, y c ~[010], 2c«[001] 

x t =1/V£[ll0], y,=l/V2[110], Z!-[001] 

x 2 =)/7l3[320], y 7 =Wl3[230], z 2 =[00i] 

The coordinate system R 2 is deduced from the coordinate 
system R x by a rotation of [001] axis through an angle 
35 y = 11.31° (see diagram 2b). The relationships between 
the microscope coordinate system and the crystal 
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coordinate systems depends on the orientation of the 
specimen in the microscope. 

Point (i) : the thin lamella with practically parallel 
5 faces was produced with an FIB (see comments 1 and 2). 
The normal to the faces was chosen to be very close to 
the yi direction (figure 2b) (but here also a different 
geometry could be chosen, for example y 2 ) . In the first 
series of experiments, the thickness t of the lamella 
10 was t = 320 nm. This thickness was measured by a 
relatively conventional convergent beam technique - see 
the book "Electron MicroDif fraction" by J.C.H. Spence 
and J.M. Zuo (Plenum Publishing Corporation) . 

15 Comments: 1) the departure from parallelism was 
measured by the convergent beam and energy loss 
techniques. These techniques are relatively 

conventional (cf. the aforementioned book "Electron 
MicroDif fraction" by J.C.H. Spence and J.M. Zuo with 

20 regards to the convergent beam and the book "Energy- 
Filtering Transmission Electron Microscopy" by Reimer, 
published by Springer Verlag in the case of energy 
loss) . The thickness of the lamella was therefore 
measured at various distances from the upper layer. It 

25 was found that the upper lamella made an angle of 1.15° 
to the lower lamella. In our feasibility demonstration, 
this angle was neglected and the lamella was considered 
to have parallel faces. However, it would have been 
possible to introduce this angle into the finite- 

30 element calculations. 

2) The FIB is also capable of producing lamella with 
better parallelism. 

35 Point (ii) : The CBED patterns were' taken along (or very 
close to) the orientation direction y 2 = [230], that is 
to say y 2 was parallel to z 0 . This is the direction used 
in the STREAM project (but other viewing directions are 
also possible. Forming CBED patterns in various 
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directions would increase the number of experimental 
data points) . Typically, the size of the electron beam 
was taken to be equal to 0.4 nm, the beam angle being 
around 15 mrad. During a scan, from 50 to 100 CBED 
5 patterns are taken every 4 nm in a direction 
perpendicular to the surface (the z 2 direction) , 
starting from the surface, but only around ten patterns 
are used in the calculations. Figure 3 indicates the 
position along the z 2 axis where five of these 
10 experimental patterns were obtained. 

Various types of CBED pattern may be defined: 

a) far from the interface between parts A and B, 
the silicon substrate is, with the measurement 

15 precision, considered as undeformed (for example, the 
-Z2 = 300 nm pattern in figure 3a) . This may serve as a 
reference for determining the experimental parameters, 
as is done in the STREAM project; 

b) by approaching the interface, the Holz lines 
20 remain fine, but start to shift slightly - the crystal 

parameters of the substrate undergo very slight 
modifications. This is the effect measured and 
quantified in the STREAM project; 

c) closer to the interface, the Holz lines broaden 
25 (for example, the -z 2 = 139 nm pattern in figure 3a) . 

We will call them Holz bands. This is the effect that 
is measured and quantified below; and 

d) even closer to the interface, the Holz lines 
become too broad and too weak: the CBED patterns no 

30 longer contain enough details to be quantified (for 
example, the -z 2 = 70 nm pattern in figure 3a) . 

Point (iii) : 

35 Once this series of measurements has been carried out 
at the thickness t = 320 nm, the lamella is thinned 
further down to a thickness t = 300 nm. A second series 
of measurements is then carried out. 
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Fourteen patterns were chosen for the first thickness 
and ten for the second (see figure 5) . This number is 
sufficient to reproduce the changes in the rotation. 

5 Point (iv) : 

The JEMS software (P. Stadelmann, CIME-EPFL CH1015- 
LAUSANNE) was used to reproduce the position of the 
Holz lines of the CBED patterns far from the interface 

10 (type a CBED patterns), but other software packages or 
the theory presented in the aforementioned book by 
Spence and Zuo could be used. If these Holz lines of 
the silicon substrate - the reference crystal - are 
transferred onto the type-c CBED patterns, these lines 

15 lie in the middle of the Holz bands (the lamella is 
homogeneous and symmetrical) . The positions of one of 
the ends of these Holz bands are simulated in the JEMS 
software by tilting the reference crystal by a rotation 
of x 2 //[3-20] axis through an angle -0 ma x (the lines in 

20 figure 5b) . The positions of the second end are 
simulated in the JEMS software by tilting the reference 
crystal, by a rotation of x 2 //[3-20] axis through an 
angle 0 max (the lines in figure 5b) . A rotation of angle 
AG = 29 raax therefore allows the variable widths A0 g of 

25 the Holz bands to be reproduced to a good 
approximation . 

Before the thinning operation by the FIB technique, the 
silicon substrate is strained, but generally only 

30 slightly deformed. After thinning, the strains produced 
by the interface between parts A and B are relaxed. At 
each coordinate point (yi,zi), the finite elements 
therefore give a displacement R(yi,Zi) of components u 
and v along the y x and Zi axes. This displacement may be 

35 decomposed' into a local translation, a pure deformation 
e(yi,Zi) and a local rotation of x x axis and through an 
angle p (yi, z 1 ) : 



The simulations in the JEMS software indicate that it 
is these local rotations P(yi,zi) along the electron 
beam (the z 0 direction) that, are the main cause of the 
broadening of the Holz lines into Holz bands. A 
rotation of Xi axis and through an angle P(yi,zi) may be 
decomposed into three rotations of x 2 axis through an 
angle 0(yi,Zi), of the y 2 axis through an angle 9'(yi,Zi) 
and of the z 2 axis through an angle 0 1 ' (y x , zi) . To a 
first order, only the first of these three rotations is 
important, giving the equation: 



In particular, it has also been shown that the angle 
©max is to a first order equal to 0.98p max , where p max is 
the maximum value of the angle P(yi,Zi) along the 
electron beam, that is to say along y 2 , or, which 
amounts to the same thing, to the maximum value of the 
angle p(yi,Zi) along the y x direction since our lamella 
is periodic in Xi : 



In this feasibility demonstration, this approximate 
property is used to optimize the constants of the 
material. Each CBED pattern taken at a depth z 2 is 
therefore characterized by a single parameter 
A9(z 2 ) = 29 max (z 2 ) (figure 3b). 

The curves shown in figure 6 plot the experimentally 
measured A9 values and the 0.98AP values (where 
AP = 2p max ) determined from the finite-element 
calculations . 



B(yi,zi) = O.980(yi,zi) . 



9 max (zi) = 0.98p max ( Zl ) = 0.98 




The calculated 
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curves in 
the NiSi 
reproduce 



figures 6a and 6b use 
material taken from the 
the experimental curves 



better, the parameters of the NiSi material are varied 
in the calculation of the calculated curves shown in 
figures 6c and 6d. By trial and error, the constants 
for the NiSi layer (Young's modulus E and Poisson's 
ratio v) are obtained (using the isotropic 
approximation) . 

A few details will now be given about the use of finite 
elements. A commercial software was used and the plane 
deformation approximation was used. 

Three particular states are important: 



state 



it is assumed that there exists a 



temperature T 0 = T x + AT at which the silicon substrate 
and the NiSi film would be coherent and unstrained 
(this temperature is not necessarily obtainable - it 
could be referred to as the "imaginary coherence 
temperature") ; 

state 2 : the thin lamella of thickness t, coherent 
at temperature T 0 , is cooled to the experimental 
temperature Ti = T 0 -AT. The strains are partly relaxed 
at the surface of the thin lamella. 

The classical laws of solid mechanics were used to 
determine the displacements R(yi,zi). In particular, the 
equation between the stresses and strains is given by: 



(l + v)(l-2v) 



I- v v 0 
v 1-v 0 
0 0 l-2v 



V 



E 



d-2v) 



a 



AT 



where AT = T 0 -T 1# E is the Young's modulus, v is the 
Poisson's ratio, and a is the thermal expansion 
coefficient. These equations apply in the Si substrate 
(E A = 156 GPa, v A = 0.277 and a A = 2.6 x 10" 6 K" 1 ) and in 
the thin NiSi film (E B , v B , a B and AT were optimized by 
minimizing the distance % between the experimental and 
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calculated curves (see figure 6)) using the assumption 
of plane strain. The strains were calculated in the 
plane defined by the direction of the incident 
electrons, yi//z 0 and the z 1 axis; and 

State 3 : before extraction of the lamella, the 
system was much more strained as the strains are not 
relaxed via the surfaces. Once the parameters (for 
example the elastic coefficient and the thermal 
expansion coefficient) have been determined, the 
strains in the integrated circuit before thinning can 
be determined by means of the present method. The 
strains in the transistor are calculated using the 
optimized constants: these are the strains taken in the 
middle of a lamella of very large thickness t 0 or those 
of a periodic lamella of infinite thickness. 

Comments : 

Ri) as regards point (iv), in this feasibility 
demonstration, the elastic constants are adjusted 
20 manually, although automatic optimization is possible. 

R 2 ) as regards point (iv), we also used the approximate 
equation: 9 max = 0.98p max , that is to say only the total 
width A0 g of the Holz bands was taken into account. A 
25 very accurate but longer adjustment would consist in 
reproducing the Holz bands in their entirety. Each CBED 
pattern would then no longer be characterized by a 
single value A0, but by a set of profiles identical to 
those shown in figure 7b. 

30 

These profiles represent in fact the function 1-A* (s) 

where A g (s) is the amplitude of the wave function 
diffracted by the planes of index g, and s measures the 
deviation from Bragg 1 s law for these planes g. As 
35 explained in the aforementioned book by Spence and Zuo, 
s and 0 are connected via a simple formula and various 
approximations may be made in order to calculate A g (s). 
It has been shown, for the fine Holz lines of CBED 
patterns taken in the y2//[230] direction, a two-wave 



10 
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approximation or a kinematic approximation gives the 
same physical result. In the kinematic approximation: 

h z 

5 The parameter to be minimized would in this case be the 
sum of the distances between curves of the type shown 
in figures 7b and 7c. The advantage of the above 
formula is that not only the rotations but also the 
pure deformations are taken into account in the R(y2,z 2 ) 
10 term (even if it has been shown that the essential 
effect in the broadening of the Holz bands is that due 
to the local rotations) . 

R 3 ) points (i) and (ii) have also been made on specimen 
15 2 described in part 3. 

R 4 ) at different temperatures, the principle of the 
method remains similar, and is relatively conventional 
in the measurement of strains. Introducing measurements 
20 at different temperatures amounts to increasing the 
number of experimental data points and the number of 
parameters to be optimized: the methodology is quite 
similar . 

25 The flowchart shown in figure 9 gives a summarizing 
view of the way in which the method of the invention 
that has just been described in detail is implemented. 

Step El of the method is the cutting of at least one 
30 lamella comprising part of the substrate and part of 
the film placed on its surface. This results in a 
relaxation of the strains (E2) and therefore a 
deformation of said lamella. At step E3, at least one 
deformation parameter (preferably a rotation P) is 
35 measured corresponding to a plurality of points on the 
lamella, at various depths from the film/substrate 
interface . 
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Independently of steps El - E3, modeling, typically by 
finite-element analysis, of the strained lamella, that 
is to say before it is cut, is carried out. To do this, 
a first estimate of the mechanical properties of the 
5 thin film, such as its Young's modulus E 1 , its 
Poisson's ratio v' and its coherence temperature T ! 0 , is 
made (step E4). Next, at step E5 the relaxation of the 
strains is modeled. This may be carried out by 
replacing imposed (zero) displacement boundary 
10 conditions with free boundary conditions. At step E6, 
the expected values s f and P f of the deformation 
parameter or parameters corresponding to the points 
where the measurements of step E3 were taken are 
determined on the basis of this simulation. 

15 

Preferably, but not necessarily, steps El - E6 are 
repeated for a plurality of lamella of different 
thicknesses, or more generally lamella having a 
different geometry, and/or at different temperatures. 

20 

The root mean square error between the measured and 
expected values is calculated (E7) and minimized by 
varying the estimates of the mechanical properties of 
the thin film (E9) and iteration of steps E4 - E7, 
25 until convergence is reached (E8). 

In this way, an optimum estimate, within the root mean 
square error, of the mechanical properties of said film 
is obtained (E10) . 

30 

Optionally, the information thus obtained is used, in 
step Ell, to calculate a deformation and/or strain 
state of the thin film and/or of the substrate. 

35 The above description has been essentially limited to 
the case of a substrate on which a thin film has been 
deposited and to measurements made on specimens formed 
by a lamella having approximately mutually parallel 
faces perpendicular to the interface between the 
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substrate and the film. However, the method also 
applies to specimens having a different shape, for 
example a wedge shape. What is important is that the 
geometry of this specimen be accurately known, in order 
5 to be able to carry out the modeling, and that it 
allows strain relaxation. 

Likewise, when measurements are carried out on a 
plurality of specimens in order to obtain a better 

10 estimate of the mechanical properties of the film, the 
mechanical property of said specimens be varied, and 
not just their thickness. For example, it is possible 
to use a series of wedge-shaped specimens having 
different angles between their faces and/or a different 

15 angle of inclination relative to the film/substrate 
interface . 

As the characterization of a surface layer was of 
particular interest, the various measurements of 

20 deformation parameters were made at different depths 
from the film/substrate interface. More generally, in 
the case of inclusions, fibers or filaments, the term 
"depth" no longer applies but simply the distance from 
the interface between the two phases in question. 

25 Furthermore, certain of the measurement points may, at 
least in some cases, lie within the nanoscale or 
microscale phase to be characterized. 

As regards the measurement technique, methods other 
30 than CBED may be used, such as for example LACBED 
(Large-Angle Convergent Beam Electron Diffraction) . 

In the example in question, the only deformation 
parameter considered was the angle of rotation p, 
35 determined by measuring the broadening A0 of the Holz 
lines. The displacement of the Holz lines, associated 
with pure deformations, could also be taken into 
account. More generally, the method of the invention 
may include the use of one or more deformation 
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parameters determined from various directly measured 
quantities . 



